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The excitation spectra in the deformed nucleus ^aorpj^ were studied by means of the (p,t) reaction, 

. . . using the Q3D spectrograph facihty at the Munich Tandem accelerator. The angular distributions 

0^ ' of tritons are measured for about 200 excitations seen in the triton spectra up to 3.3 MeV. Firm 

I 0"*" assignments are made for 16 excited states by comparison of experimental angular distributions 

■ with the calculated ones using the CHUCK code. Additional assignments are possible: relatively 

(N . firm for 4 states and tentative also for 4 states. Assignments up to spin 6 are made for other states. 

I Sequences of the states are selected which can be treated as rotational bands and as multiplets of 

^ , excitations. Experimental data are compared with interacting boson model (IBM) and quasiparticle- 

' phonon model (QPM) calculations. 
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^ ■ I. INTRODUCTION 

A full microscopic description of low-lying excitations in deformed nuclei has eluded theoretical studies to date. 
Along with the interplay of collective and single-particle excitations, which takes place in deformed rare earth nuclei, 
^ ■ additional problems arise in the actinide region because of the reflection asymmetry [1]. Evidently the nature of the 
, first excited 0+ states in the actinide nuclei is different from that in the rare earth region where they are due to the 
' quadrupole vibration [2]. Octupole degrees of freedom have to be important in the actinides. One has then to expect a 
complicated picture at higher excitations: residual interactions could mix the one-phonon and multiphonon vibrations 
of quadrupole and octupole character with each other and with quasiparticle excitations. Detailed experimental 
information on the properties of such excitations is needed for comparison with theory. On the experimental side, the 
(p,t) reaction is very useful. On the theoretical side, a microscopic approach such as the quasiparticle-phonon model 
(QPM) is necessary, in order to account for the number of states detected and to make detailed predictions on their 
properties . 

Our previous paper [3] concentrated on the excited 0+ states in the actinide nuclei ^^^Th, ^^°Th and ^32^ studied in 
the (p,t) reaction. This interest in the 0+ excitations in the deformed actinide nuclei was stimulated by the observation 
of multiple L = transfers in the (p,t) reaction leading to the excited states in the odd nucleus ^^^Pa [4] as well 
Ci ^ as to the excited states in the medium heavy even nuclei ^"^^Nd [5], ^'^^Sm [6], ^^^'^■^*Ba [7]. After comprehensive 
information on 0+ excitations was obtained in the actinide nuclei [3] , systematic studies of the same type were made 
on many nuclei between ^^^Gd and ^^^Hg [8] and revealed a large number of excited 0"*" states, whose nature is not 
yet understood. The (p,t) reaction, however, gives much more extensive information on specific excitations in these 
nuclei, which was not analyzed previously [3]. An attempt to obtain such information was made for ^^^Er after using 
a high-resolution experimental study with the (p,t) reaction [9]. In this paper we present the results of a careful 
and detailed analysis of the experimental data from the high-precision, high-resolution study of the ^'^^Th(p,t)^'^°Th 
reaction carried out to get deeper insight into excitations in ^^°Th including the nature of the 0"*" excitations. 

Experimentally accessible thorium nuclei span a wide and interesting range of isotopes: from the semi-magic ^^^Th 
to quadrupole-deformed ■^■^^Th. An especially interesting region is that around A=228, where the even isotopes 
^^^Th, ^^^Th and ^^°Th are considered to be octupole-deformed, octupole soft and vibrational-like, respectively [10]. 
Two of them are accessible for study by the (p,t) reaction. Besides 0+ excitations, whose number is increased in 
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FIG. 1: Triton energy spectrum from the ^^^Tli(p,t)^^°Th reaction (Ep=25 MeV) in logarithmic scale for a detection angle of 
7.5°. Some strong lines are labeled with the corresponding level energies in keV and by the spins assigned from the DWBA fit. 



comparison with the prehminary analysis in publication [3], the spins for many other states are also assigned. The 
results concerning 0+, 2+ and 4+ excitations are mainly discussed. The 6+ and negative parity states are excited 
weakly in the (p,t) reaction; therefore information on these states is only fragmentary. The results of a similar analysis 
for ^^*Th in comparison with ^^^Pa as well as for 232jj ^-^^l be presented in forthcoming papers [11]. 

In Sec. II A of this paper the details of the experimental techniques and the experimental results are given. 

After description of the DWBA analysis a preliminary treatment of the nature of some excitations is presented. In 
Sec. Ill the interpretation of 0^ excitations is discussed. The experimental data obtained in this study are compared 
with the results of calculations in microscopic approach of the quasiparticle-phonon model (QPM). The validity of 
the phenomenological approach of the interacting boson model (IBM) is also tested to account for at least the main 
features of the positive parity excitations in ^^°Th, in addition to 0+ states tested in [3]. 
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FIG. 2: Comparison of the spectra for the ^^^Th(p,t)^^°Th reaction (Ep=25 MeV) for detection angles of 12.5° (thin Une) and 
26° (thick hne) for the high energy range. The Unes corresponding to possible 0"*" excitations are labeled by their energies. 



II. EXPERIMENT, ANALYSIS AND EXPERIMENTAL DATA 
A. Details of the experiment 

A target of 100 /xg/cm^ ^^^Th evaporated on a 22 ng/cm^ thick carbon backing was bombarded with 25 MeV 
protons of an intensity of 1-2 jiA from the Tandem accelerator of the Maier-Lcibnitz-Labor of the Ludwig-Maximilians- 
Universitat and Technische Universitat Miinchen. The isotopic purity of the target was about 99%. The tritons were 
analyzed with the Q3D magnetic spectrograph and then detected in a focal plane detector. The focal plane detector 
is a multiwire proportional chamber with readout of a cathode foil structure for position determination and dE/E 
particle identification [12, 13]. The acceptance of the spectrograph was 11 msr, except for the most forward angle of 5° 
with an acceptance of 6 msr. The resulting triton spectra have a resolution of 4 7 kcV (FWHM) and are background 
free for all angles but 5° for which background from light contaminations in the region of 1100-1150 keV complicated 
the analysis for some levels. The angular distributions of the cross sections were obtained from the triton spectra at 
ten laboratory angles from 5° to 45°. 

A triton energy spectrum measured at a detection angle of 7.5° is shown in Fig. 1, which demonstrates the domina- 
tion of 0+ excitations having large cross sections at this angle. The analysis of the triton spectra was performed with 
the program GASPAN [14]. Measurements were carried out with two magnetic settings: one for excitations up to 1.6 
MeV and another for the energy region from 1.5 MeV to 3.3 MeV. For the calibration of the energy scale the triton 
spectra from the reaction ^*"'W(p,t)^^^W, ^^^W(p,t)^^'*W and ^■''*U(p,t)^'^^U were measured at the same magnetic 
settings. Some known levels in ^^°Th were included in the calibration, leading to small corrections from to 0.5 keV 
for excitations between 1.5 MeV and 3.3 MeV. In the course of our measurements we found that the Q values for the 
(p.t) reactions on the ^'^^Th, ^sor^j^ g^j^j ^-^^U targets or on the -'^^^W and ^^^W targets are in disagreement with the 
ones calculated from the Atomic Mass Evaluation [15] depending which Q value is used as the reference value. They 
are given in Table I. Notation "Ref." (the reference value) in this table means that the Q value for this nucleus as 
determined from the data in the AME was taken as a starting point in the calculations for other nuclei. If the Q 
value for the ^^^Th(p,t) reaction as derived from [15] is taken as the reference value, then the Q values for ^^°Th(p,t) 
and ^•^^U(p,t) reactions also agree with those calculated from the AME table [15]. At the same time the Q values 
for the ^^^W(p,t) and ^*''W(p,t) reactions differ considerably from the calculated one and vice versa, if the Q value 
for the ^^^W(p,t) reactions is taken as the reference value. It is not clear from these data what is the reason for this 
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TABLE I: Q values for the (p,t) reactions as determined from a comparison of the experimental triton spectra and as derived 
from the AME 2003 compilation [15]. 



Reaction 


Hei.: " in(p,t) 


Hei.: W(p,tj 


AMbi 2003 [15J 




[keV] 


[keV] 


[keV] 


232Th(p,t) 


-3076.5(27) 


-3064.2(20) 


-3076.5(27) 


23r'Th{p,t) 


-3568.5(30) 


-3556.6(20) 


-3569.0(28) 


234u(p,t) 


-4125.2(30) 


-4113.0(20) 


-4124.2(29) 


i8*5W(p,t) 


-4478.1(30) 


-4460.5(20) 


-4463.0(19) 


"*W(p,t) 


-5133.0(30) 


-5120.6(12) 


-5120.6(12) 



discrepancy. 

About 200 levels have been identified in the spectra for all ten angles and are listed in Table II. The energies 
and spins of the levels as derived from this study are compared to known energies and spins from the published data 
[16-18]. They arc given in the first four columns. The ratios of cross sections at angles 5° and 26° to that at angle 16°, 
given in the next two columns, help to highlight the 0+ excitations (large values). The column crjjj^gg gives the cross 

section integrated in the region from 5° to 45° and the column crexp./^calc S^^^^ the ratio of the integrated cross 
sections, from experimental values and calculations in the DWBA approximation (see Sec. II B). The last column 
gives the notations of the schemes used in the DWBA calculations: sw.jj means one-step direct transfer of the (j)^ 
neutrons in the (p,t) reaction; notations of the multi-way transfers used in the DWBA calculations are displayed in 
Fig. 6 (see Sec. II B). 

In order to get an indication of the possible 0+ excitations at higher energies, the triton spectra were measured for 
the angles of 12.5° and 26.0° and in the energy range from 2.8 MeV to 4.5 MeV. A precise calibration of this region 
of excitations is not possible with the spectra measured for the tungsten and uranium targets or any other targets. 
Therefore the calibration makes use of that employed for the energy region 1550 - 3300 keV with a linear shift of the 
energy and channel number determined from the position of the known levels observed in both energy regions. The 
measured spectra are compared in Fig. 2. The ratio of the line intensities at 2808 keV in two spectra can serve as 
a reference for the 0+ state assignments. Peaks that may correspond to the excitation of 0+ states are marked by 
dashed lines. Unfortunately, close values of the excitation ratios are valid also for 6+ and 3+ excitations as shown 
by the states at 3052 and 3072 keV. Nevertheless, at least 10 candidates for 0+ excitations are revealed in the region 
from 3 to 4 MeV. Measurements at the angle of 5° could confirm or reject these assumptions. 

TABLE 11: Energies of levels in ^'^"Th, the level spin assignments 
from the CHUCK analysis, the (p,t) cross sections integrated over 
the measured values and the reference to the schemes used in the 
DWBA calculations (see text for more detailed explanations). 



Level energy [keV] 
This work [16-18] 


r 

[16-18] This work 


Cross section ratios 
(5716°) (26716°) 


'^integ. 

[nh] 


Ratio 
""expt./'^'calc. 


Way of 
fitting 


0.1 2 


0.0 


0+ 


0+ 


8.02 


5.51 


195.68 


10.7 


sw.gg 


53.2 2 


53.20 2 


2+ 


2+ 


1.57 


0.27 


52.53 


11.0 


mla 


174.0 2 


174.10 3 


4+ 


4+ 


0.87 


0.50 


9.94 


2.6 


ml a 


356.3 2 


356.6 5 


6+ 


6+ 


0.40 


0.75 


7.39 


2.8 


m2d 


508.0 3 


508.15 3 


1" 


1" 


0.63 


0.00 


0.87 


0.6 


m2a 


571.7 2 


571.73 3 


3" 


3" 


0.33 


0.44 


4.04 


0.8 


m3a 


593.8 3 


594.1 5 


8+ 


8+ 


1.05 


0.00 


0.37 


0.4 


m2c 


635.1 2 


634.88 5 


0+ 


0+ 


26.91 


8.92 


47.79 


250 


sw.ii 


677.6 2 


677.54 5 


2+ 


2+ 


0.85 


0.55 


25.13 


3.5 


mla 


686.0 10 


686.7 


5" 








< 0.2 






775.2 4 


775.5 3 


4+ 


4+ 


0.30 


0.53 


7.15 


1.3 


mla 


781.4 2 


781.35 3 


2+ 


2+ 


0.62 


0.46 


69.07 


6.2 


sw.gg 


825.6 3 


825.66 5 


3+ 


3+ 


0.14 


0.79 


1.10 


2.5 


m2a 


852.7 4 


851.88 3 


7" 








0.42 






884.2 4 


883.9 2 


4+ 


4+ 


0.51 


1.27 


4.09 


0.8 


mla 



TABLE II: Continuation 



Level energy [keV] 


r 


Cross section ratios 


''^integ. 


Ratio 


Way of 


This work [16-18] 


[16-18] 


This work 


(5'V16°) 


(26°/16°) 


[/ib] , 


'^cxpt./'^cak:. 


fitting 


923.3 5 


923.0 1! 


6+ 


6+ 


1.66 


0.96 


u.ou 


u.z 


Trt 9r1 


952.6 5 


951.88 5 
955.1 2 


1" 
5+ 


1", 
or (0+) 


1.17 


1.76 


n 71 

U. ( i 


u.uo 
u.u 


ml 

in la 
sw.ij^ 


972.1 5 


971.72 5 


2" 


2" 


0.22 


0.69 


U.O 1 


o.u 


m9Q 

mza 


1011.6 5 


1009.58 5 


2+ 


2+ 


0.54 


0.34 


14.34 


2.1 


sw.gg 




1012.51 5 


3" 


3" 








6.0 


sw.gg 


1040.0 7 


1039.6 2 


6+ 


6+ 


0.80 


1.07 


n Q/1 


U.O 


m9a 

inza 


1052.0 7 


1052.31 5 


3+ 


3+ 


0.00 


3.33 


0.34 


0.8 


m2a 


1065.9 8 






4" 


0.99 


0.69 


0.30 


0.45 


m3b 


1079.4 8 


1079.21 3 


2" 


2" 


0.00 


0.00 


0.09 


0.1 


m3b 


1108.7 5 


1107.5 5 
1109.0 i 
1117.5 3 


4+ 
5~ 
8+ 


4+ 


1.20 


0.87 


O.Oo 


1 9 


m9Q 

mza 


1125.6 5 


1127.76 5 
1134.4 2, 


3" 
7+ 


(1"). 
or (0+) 


4.54 


1.06 


2.10 


0.1 

i.O 


mla" 
sw.jj" 


1148.0 9 


1176.1 5 


5+ 














1184.8 S 


1196.8 
















1241.2 9 


1243.3 3 


8+ 








<0.2 






1256.0 9 


1255.5 5 


6+ 








<r'n 9 






1259.2 6 






(3-) 


0.97 


0.34 


U.04I: 


U.O 


sw.gg 


1283.6 6 






(5-) 


1.08 


0.00 


U.Zo 


u.u ( 


moD 


1297.8 6 


1297.1 J 


0+ 


0+ 


4.37 


2.53 




i.O 


sw.ig 


1322.3 5 






(3-) 


1.23 


0.54 


1 IT 


i.D 


sw.gg 


1337.2 5 


1349.3 4. 


7+ 


4+ 


1.37 


1.14 




^ 9 


sw.gg 


1359.5 7 






(2+) 


0.00 


0.41 


n AS 

U.4I:0 


OPi 

u.uo 


sw.gg 


1376.6 7 


1375.3 1 


1,2+ 


l+,5- 


0.00 


0.00 


n 9/1 


n 9 /n 1 
u.z/ U. i 


moD 


1401.5 5 


1400.9 1 


2+ 


2+ 


0.49 


0.38 


1 s 
o. lo 


U.O 


sw.gg 


1420.4 5 






(3+) 


0.33 


0.70 


A 7/1 


i.U 


m 9r, 

mza 


1440.4 8 






(3+) 






9 






1447.9 5 






0+ 


10.37 


6.86 


1 70 


0^ 
u.uo 


sw.gg 


1485.6 5 


1485.6 1 




4+ 


1.44 


0.93 


4.59 


0.6 




1496.0 10 












<0.2 






1507.4 5 






4+ 


1.29 


2.28 


n "17 


07 


m2a 


1524.8 5 






2+ 


0.66 


0.00 


0.53 


0.05 


bW.gg 


1566.2 6 








1.80 


0.29 


0.31 


0.15 


m2b 


1574.5 6 


1573.5 a 


l(-),2+ 


(2-) 


0.34 


0.70 


2.14 


o.o 


m2a 


1584.7 6 






(4-, 5+) 


1.77 


1.90 




/1 9 7*^^ 
4Z/ 00 


niza 


ioyU.z 


iOoy.o Z 


0+ 


0+ 


8.48 


6.04 


11. DO 


U.Zo 


sw.gg 


1594.7 « 








1.11 


0.58 


9 99 


-L.O 


mou 


1601.2 






(3-) 


1.46 


1.14 


n 74 


S 
U.O 


sw.gg 


1612.1 10 






(4-, 5+) 


1.26 


2.25 


u.ou 


O.o/ 1 


m9n 


1618.7 9 






(4-, 5+) 


1.44 


2.09 


n 97 

U.Z ( 


/O Q 
o.u/ u. y 


m2a 


1630.1 7 


1628 a 




2+ 


1.43 


0.45 


4.49 


0.6 


mla 


1639.3 6 


1638.5 2 


(2,0+) 


0+ 


5.68 


4.64 


8.32 


0.2 


sw.gg 


1653.2 






(6+) 


0.96 


1.75 


0.88 


0.004 


m3a 


1668.2 7 






4+ 


1.54 


1.85 


1.91 


0.3 


sw.gg 


1679.1 7 






2+ 


1.72 


0.36 


3.36 


0.45 


mla 


1683.3 7 






(4-) 


0.19 


1.95 


2.20 


34 


m2a 


1694.9 7 


1695.7 1 


l(-),2+ 


(4+) 


1.09 


1.41 


3.57 


0.35 


sw.gg 



TABLE II: Continuation 



Level energy [keV] 


r 


Cross section ratios 


''^integ. 


Ratio 


Way of 


This work [16-18] 


[16-18] 


This work 


(5716°) 


(26°/16°) 


[/i.b] . 


^cxpt./^calc. 


fitting 


1708.8 


8 




2+ 


0.22 


0.18 


U.oo 


U.Uo 


sw.gg 


1723.5 


7 




(4+) 


0.55 


1.35 


2.31 


0.3 




1745.3 


1744.9 1 




0+ 


2.80 


2.59 


1.18 


0.10 


SW.Jg 


1 'yen T 


o 

O 




(3-) 


1.76 


0.64 


QQ 


u.o 


sw.gg 


1 Tan o 


o 
O 




(4+) 


0.73 


0.43 




n 9 


mla 


i ( D9.D 


O I (({).( 1 


1,2+ 


(4+) 


0.52 


0.69 


u.oy 


U. J-tJ 


sw.gg 


1 77/1 1 


Q 1 77^ 9 1 


1,2+ 








<0.2 








1789.4 5 


l(-),2+ 














1793.1 


6 




(5-) 


1.19 


1.19 


2.43 


3.5 




18UZ.0 







0+ 


8.04 


5.38 


2.38 


0.05 






1810.7 1 


1,2+ 














1812.0 


8 




4+ 


1.71 


0.96 


1.26 


0.13 




1824.9 


7 




(6+) 


0.45 


0.96 


1.15 


0.15 




1840.0 


8 1839.6 2 


l(-),2+ 


2+ 


1.74 


0.70 


2.37 


0.25 




1851.4 


7 1849.6 1 


(2+) 


(3-) 


1.07 


0.93 


1.73 


1.5 




1859.3 


7 1858.2 6 




(3-) 


1.41 


1.05 


1.74 


1.6 


QiTj O'er 


1868.9 


7 




(0+) 


1.99 


1.60 




U.UO 


sw.gg 








+(6+) 








n 9fi 

U.ZD 


mza 


1887.0 


9 




(2+) 


0.84 


0.68 


n 7(\ 

U. 1 D 


U.UD 


sw.gg 




1902.7 1 


1,2+ 














1910.0 


9 




(6+) 


0.51 


1.13 


1 8/1 


n 9t^ 


sw.gg 


1914.7 


9 




(1") 


0.91 


0.57 


lis 

1. lo 


n 8 
u.o 


m9a 


1926.0 


7 




4+ 


2.18 


0.90 


1 71 


u.o 


m9a 


1931.1 


8 






0.81 


0.50 




u.o 


m9Q 


1939.8 


11 




(1^,1+) 


0.36 


0.68 


u.yo 


u.o 


m9Q 
IIlZd> 


1947.0 


6 




4+ 


1.18 


0.89 


O.Oi 


u.^ 


mla 




1949.8 1 


1,2+ 














1956.4 


6 




2+ 


0.39 


0.47 




u.oo 


sw.gg 


1967.1 


7 1966.9 2 


1,2+ 


2+ 


0.27 


0.66 


^ 7(\ 
0. 1 D 


n 1 8 

U. io 


sw.gg 


1972.0 


9 1973.4 2 


(l+,2+) 


2+ 


0.58 


0.39 




U.UO 


sw.gg 


1985.4 


8 




(5-) 


0.00 


0.52 


n t^i 


n 8 
u.o 


sw.gg 


2001.6 


8 2000.9 1 


1,2+ 


(3-) 


1.28 


0.83 


1 1 

1.10 


1 1 


sw.gg 


2010.3 


6 2010.1 2 


1,2+ 


2+ 


0.26 


0.38 


o.oy 


u.o 


sw.gg 


2017.3 


7 




(3-) 


1.03 


0.73 


1 7^ 
1. 1 


i.O 


sw.gg 


2025.6 


6 2024.7 2 


l+,2+ 


2+ 


0.28 


0.39 




n 9 

U.Z 


sw.gg 


2032.8 


7 




4+ 


1.07 


1.13 


O.Dy 


u.o 


sw.gg 


2039.1 


7 




4+ 


2.24 


1.85 


zl f^l 


n zl 

U.-rt 


sw.gg 


2048.7 


7 




(4+) 


1.00 


1.19 


1 PiQ 


n 1 
u. 10 


sw.gg 


2060.9 


12 




(3-) 


2.05 


0.55 


U.o ( 


n n9 
u.uz 


ni3a 


2073.2 


8 




(8+) 


0.00 


1.82 


1 /in 


n /I 


sw.gg 


2074.9 


8 




(4+) 


0.85 


0.51 


1 1 

1. ID 


n 9 
u.z 


mla 


2085.9 


8 




(4+) 


1.75 


1.26 


1 (^n 

l.DU 


n 1 
U.lo 


sw.gg 


2093.9 


7 




0+ 


4.44 


2.17 


D.yo 


1 Q 
lo 


sw.ii 


2102.0 


7 




4+ 


2.30 


1.10 


i.zD 


u.lo 


sw.gg 


2118.4 


6 




4+ 


2.32 


1.14 


7 7n 


u.yu 


sw.gg 




2122.8 1 


1,2+ 














2130.7 


7 




2+ 


0.51 


0.76 


O.ZO 


n 9 
u.z 


sw.gg 




2133.2 2 
















2137.9 


7 




2+ 
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FIG. 3: CHUCKS one-step DWBA calculations of the angular distributions for different (j)'^ transfer configurations. The lines 
are marked with the (j)^ transfer configuration and a scaling factor introduced to allow comparison with the {2gQ^2)^ transfer 
configuration. 
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For levels at 952.6 and 1125.6 keV the angular distributions can be fitted by the DWBA calculations as spin 1 or spin 0+. 
For the levels at 1868.9 and 2422.7 keV the angular distributions are fitted by the DWBA calculations as doublet lines. 



B. DWBA analysis 

The spins of the excited states in the final nucleus ^'^'^Th were assigned via an analysis of the angular distributions 
of tritons from the (p,t) reaction. In a previous publication [3] the angular distributions for 0+ excitations were 
demonstrated to have a steeply rising cross section at very small reaction angles and a sharp minimum at a detection 
angle of about 14°. The angular distribution for the 0+ ground state of ^^°Th has such a shape. This pronounced 
feature helped us to identify these states in complicated and dense spectra without fitting experimental angular dis- 
tributions. No complication of the angular distributions was expected, since the excitation of 0"*" states predominantly 
is a one-step process. This is not the case for the excitation of states with other spins, where multi-step processes 
could play a very important role. 

The identification of other states is possible by fitting the experimental angular distributions with those calculated 
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TABLE III: Optical potential parameters used in the DWBA calculations. 
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"According to [21]. 
'According to [22]. 



in tlic distorted- wave Born approximation (DWBA). A problem arising in such calculations is that we have no prior 
knowledge of the microscopic structure of these states. We can assume, however, that the overall shape of the angular 
distribution of the cross section is rather independent of the specific structure of the individual states, since the 
wave function of the outgoing tritons is restricted to the nuclear exterior and therefore to the tails of the triton form 
factors. To verify this assumption DWBA calculations of angular distributions for different (j)^ transfer configurations 
to states with different spins were carried out. The coupled channel approximation (CHUCK3 code of Kunz [19]) was 
used in these calculations. 

The shape of the calculated angular distributions depends strongly on the chosen potential parameters. We used 
parameters suggested by Becchetti and Greenlees [20] for protons and by Flynn et al. [21] for tritons. These parameters 
have been tested via their description of angular distributions for the ground states of ^^^Th, ^^^Th and ^32^ [3]. They 
are listed in Table III. Minor changes of the parameters for tritons were needed only for some 3~ states, particularly 
for the state at 571.7 keV. For these states the triton potential parameters suggested by Becchetti and Greenlees [22] 
were used (the last column in Table III). For each state the binding energies of the two neutrons are calculated to 
match the outgoing triton energies. The corrections to the reaction energy are introduced depending on the excitation 
energy. The best reproduction of the angular distribution for the ground state was obtained for the transfer of the 
(259/2)^ configuration in the one-step process. This orbital is close to the Fermi surface and was considered as the 
most probable in the transfer process. Other transfer configurations that might be of importance are (liii/2)^ and 
(1^15/2)^: since these orbitals are also near the Fermi surface. All other configurations for the orbitals in the vicinity 
of the Fermi surface were tested in the one-step calculations. As one can see in Fig. 3, the shape of the angular 
distributions depends to some degree on the transfer configuration, the most pronounced being found for the 0+ 
states. However, the main features of the angular distribution shapes for 2'^ and 4"^ states are not dependent on the 
transfer configurations. Therefore the (2^19/2)^ configuration was used in the calculations for the majority of excited 
states. For some 0"*" excitations the configurations (lHi/2)^ and (1^15/2)^ alone or in combination with (2^9/2)^ were 
needed to obtain better reproduction of the experimental angular distributions. 

Results of fitting the angular distributions for the states assigned as 0"*" excitations are shown in Fig. 4. The 
angular distributions in the two first columns measured with good statistical accuracy are believed to give firm 0"*" 
assignments. Three possible transfer configurations (2^9/2)^, (liii/2)^ and (lji5/2)^ have been used to get the best 
fit to the experimental data. These configurations are listed in the last column in Table II. The assignments to the 
states of 1745.3, 2368.5 2426.4 and 2596.4 keV are considered as relatively firm, because the shape of the angular 
distribution is fitted perfectly by the calculations although there is a limited statistical accuracy of the experimental 
data. The angular distribution for the transition at the excitation energy 1868.9 keV is fitted as a sum of two angular 
distributions for transfer to a 0+ and a 6+ state (doublet line). A similar situation is assumed for the energy 2422.7 
keV, i.e. a superposition of two angular distributions for 0+ and 4+ states. In both cases assignments are tentative. 

The angular distributions for the energies 952.6 and 1125.6 keV are a special case. A state with spin 1~ is known at 
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angle [deg] 

FIG. 4: Angular distributions of assigned O""" states in '^^°Th and their fit with CHUCK3 one-step calculations. The {ij) transfer 
configurations used in the calculations for the best fit are given in Table II. The first two columns on the left correspond to 
firm assignments and the column on the right to tentative assignments. 



an energy 951.9 keV, however the angular distribution of the tritons at an excitation energy of 952.6 keV is completely 
different from that for a known 1~ state at 508.1 kcV. Nevertheless, it can by fitted satisfactorily by an inclusion 
of one-step and two-step excitations. At the same time this angular distribution can also be fitted by a calculation 
for a 0^ excitation mainly by the (lji5/2)^ transfer configuration with a small admixture of the (liii/2)^ transfer. 
A spin of 0+ is assigned to the state with closely lying energy of 927.3 keV in ^'^^U, studied in the alpha decay of 
236py j23]. Therefore, we have not excluded spin assignment 0+, though the present information is not sufficient to 
solve the problem. A similar angular distribution is observed for the excitation at 1125.6 keV. But in this case the 
known level at the closely lying energy of 1127.8 keV has spin of 3~. The CHUCK calculations give a satisfactory fit 
for excitations of the 1~ and 0+ states but not for 3~ state. It is interesting to note that the IBM calculations (see 
below) predict 0+ excitations at close energies for both states. 
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FIG. 5: Experimental distribution of the (p,t) strength integrated in the angle region 0° - 45° for 0+, 2+ and 4+ states in ^^^Th. 
The levels identified reliably are indicated by filled circles and those identified tentatively are indicated by open diamonds. 
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FIG. 6: Schemes of the CHUCKS multi-step calculations tested with spin assignments of excited states in ^^''Th (see Table II). 



Thus we can make firm spin assignments for 16 states, relatively firm assignments for 4 states and tentative 
assignments for 4 states, in comparison with 14 states found in the preliminary analysis of the experimental data [3]. 

It is difficult to estimate the additional number of 0+ excitations above 3 MeV from the spectra at 12.5° and 26° 
in Fig. 2. The intensities of the lines corresponding to 0^ states have to be very low for 12.5° and much higher for 
26°. The lines fulfilling this condition are labelled by their energies. Unfortunately, a similar condition is fulfilled 
also for 6+ excitations as well as for less probable 2~, 3+ states. To be sure of the presence of 0+ states in this high 
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FIG. 7: Angular distributions of some excited states of natural parity and their fit by the CHUCKS calculations. The (y) 
transfer configurations and schemes used in the calculations for the best fit are given in Table II. 



energy region at least a spectrum for the most forward angle has to be measured. But there is an impression from 
the measured spectra that the density of 0+ states decreases for energies above 3 MeV (or else that the cross section 
of such excitations is very low and they are hidden in very dense and complicated spectra). The maximum density 
of 0^ states is observed in the interval between 2 and 3 MeV, as one can see in Fig. 5. The neutron pairing gap for 
230rpj^ is about 750 kcV and this leads to two-quasiparticle excitations around 1.5 MeV. Thus the maximum density 
of 0+ states above this energy may be caused by inclusion of such excitations, as predicted by the calculations [24]. 

Similar to 0+ excitations, the one-step transfer calculations give a satisfactory fit of angular distributions for about 
70% of the states with spins different from 0+ but about 30% of these states need the inclusion of multi-step excitations. 
Multi-step excitations have to be included to fit the angular distributions already for the 2+, 4+ and 6+ states of the 
ground state band. Fig. 6 shows the schemes of the multi-step excitations tested for every state in those cases where 
one-step transfer does not provide a successful fit. Fig. 7 demonstrates the quality of the fit of some different-shaped 
angular distributions for the excitation of states with spins higher than 0"*" by calculations assuming one-step and 
one-step plus two-step excitations. The fits for the ground state band are included in Fig. 7. Whereas natural parity 
states can be populated by one-step or one-step plus two-step mechanisms, the states of unnatural parity can be 
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FIG. 8: Angular distributions and their fit by the CHUCKS calculations for some excited states of unnatural parity, and for 
the states where our assignments of spins are in contradiction with the compilation [16]. The (ij) transfer configurations and 
schemes used in the calculations for the best fit are given in Table II. 



excited only by two-step excitations. The angular distributions and their fits for such excitations are presented in 
Fig. 8. 

The assignments of the spins resulting from such fits are presented in Table II together with other experimental 
data. Detailed explanations of this table are given in Sec. II A. Special comments are needed for the column displaying 
the ratio <Je.xpl <^cai- Since wc have no a priori knowledge of microscopic structure of the excited states and thus we 
do not know the relative contributions of the specific (j)^ transfer configurations to each of these states, these ratios 
cannot be considered as spectroscopic factors. Nevertheless, a very large ratio, such as in the case of the (I«ii/2)^ 
transfer configurations used in the calculation for the 0+ state at 635 keV, is unexpected. 

Some additional comments on Table II are needed. There are some contradictions in the assignment of the energy 
of the second 4+ level. It has been proposed from Coulomb excitation to be at 772.1 keV in [25] and at 769.6 keV in 
[26] and from the (d,pn7) reaction at 775.5 kcV [18]. The value 769.6 kcV is accepted in the compilation [16]. We see 
a weak peak at 775.2 keV on the tail of the very strong peak at 781.4 kcV as confirmation of the 775.5 kcV assignment 
[18]. There is not even a hint of a peak at 769.6 keV. Again the 8+ and 10+ levels proposed in [26] at 1251.4 kcV and 
1520.4 keV (both accepted in the compilation [16]) are in contradiction with the assignment of the 8+ level at 1243.2 
keV and with the calculated energy of the 10"*" level at 1487 keV [18]. The (p,t) study confirms the assignment for 
the 8+ level at 1243.2 keV by the observation of a weak peak at 1241.2 kcV and no peak in the vicinity of 1251.4 keV. 
After this confirmation the smooth change of the inertial parameter in the band prefers the energy 1487 keV for the 
10+ level and rejects the energy of 1520.4 keV. 

There are several levels in ^^°Th for which spins l(-), 2+ or 1,2+ were assigned from the /3 -decay and/or inelastic 
scattering and for which the (p,t) angular distributions are measured: 1573.5, 1695.7, 1770.7, 1839.6, 1849.6, 1966.9, 
1973.4, 2000.9, 2010.1 and 2024.7 keV (energies are from the compilation [16]). Angular distributions from the (p,t) 
reaction confirm spin 2+ for some of them: 1840.0, 1967.1, 1972.0, 2010.3 and 2025.6 keV (energies as determined 
from the (p,t) reaction). The (p,t) angular distributions for other states cannot be fitted for the spins given in [16]. If 
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FIG. 9: Collective bands based on the 0+, 2+, 4+, 1~, 2" and 3" excited states in '^^"Th. as assigned from the DWBA fit of 
the angular distributions from the (p,t) reaction. 

the energies 1574.5(6) keV from the (p,t) reaction and 1573.5(2) keV from [16] correspond to the same level, then a 
different assignment is suggested by the (p,t) angular distribution to this state as 2~ or 3+. In any case a 2+ excitation 
would manifest itself in the (p,t) angular distribution and therefore has to be excluded. For the states at 1694.9 and 
1769.6 keV we observe angular distributions which can be fitted by calculations for 4+ states (or even for a 0+ state 
plus a constant of about 2 /zb for the state of 1694.9 keV). Either way, the assignment 2+ can be excluded. Finally, 
the (p,t) angular distributions for the states at 1851.4 and 2001.6 keV prefer an assignment of 3~ in contradiction to 
(2+) and 1,2+, respectively, as given in [16]. 

C. Collective bands in ^^^Th 

After the assignment of spins to all excited states the sequences of states can be distinguished which show the 
characteristics of a rotational band structure. An identification of the states attributed to rotational bands was made 
on the following conditions: 

a) the angular distribution for a state as band member candidate is fitted by the DWBA calculations for the spin 
necessary to put this state in the band; 

b) the transfer cross section in the (p,t) reaction to the states in the potential band has to decrease with increasing 
spin; 

c) the energies of the states in the band can be fitted approximately by the expression for a rotational band 
E = Eq + AI{I + 1) with a small and smooth variation of the inertial parameter A. Collective bands identified in 
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TABLE IV: The sequences of states qualifying as candidates for rotational bands (from the CHUCK fit, the (p,t) 

cross sections and the incrtial parameters). More acciu-ate vahies of energies are taken from the first two columns of 
Table II. The energies taken in brackets correspond to the sequences assigned tentatively. 
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such a way are shown in Fig. 9 and are listed in Table IV (for a calculation of the moments of inertia). The procedure 
can be justified in that some sequences meeting the above criteria are already known from gamma-ray spectroscopy 
to be rotational bands, so other similar sequences are very probably rotational bands too. The straight lines in Fig. 9 
strengthen the argument for these assignments. For example the mean deviation of the experimental energies from 
the calculated rotational values for the longest newly assigned band based on the state 0+ at 1589 keV is only 1.0 
keV; for the band based on the O"*" state at 2093.9 keV it is 1.3 keV; and for the band above the 0+ state at 2268.9 
keV it is 3.9 keV. Even for the band above the 0+ state at 2426.4 keV assigned tentatively the deviation is less than 
1 keV. The observed deviations are all consistent with the stretching effect typical for rotational bands. Nevertheless 
additional information (on E2 transitions at least) is needed to confirm these assignments. 

It is worth mentioning that the tentative assignment 0+ for the state at 2426.4 keV is supported by a sequence of 3 
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TABLE V: Moments of inertia for the bands in '^^°Th as assigned from the angular distributions from the ^^^Th(p,t)^^°Th 
reaction. Results derived from the sequences having only 2 levels or assigned tentatively are given in lower part of the table. 
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other states (2+, 4+, 6+) on top of it. Three other tentative 0+ states (1745.3, 1868.9, and 2422.7 keV) have only one 
tentative state (2+) on top of them and the band sequence is not based on 7 - ray transition but on energy arguments. 
In Table V we present moments of inertia (Mol) obtained by fitting the level energies of the bands displayed in 
Fig. 9 by the expression E = Eq + AI{I +1). In upper part of the table those sequences are presented which are 
connected by known 7 - ray transitions or have at least 3 levels and in lower part of the table the sequences having 
only 2 levels or tentatively assigned are presented. The obtained Mol cover a broad range, from ~50 MeV^^ to 
~100 MeV^^. The negative parity bands based on the states with spin 1~ interpreted as the octupole-vibrational 
bands [18] have high Mol (the 1~ band at 951 keV has the largest). The 0"*" band at 1297 keV, considered as (3 
- vibrational band, has the smallest Mol. At this stage, it is difficult to make a complete correlation between the 
intrinsic structure of the bands and the magnitude of their Mol. Nevertheless one can assume also for the 0+ bands 
that the largest Mol could be related to the octupole phonon structure and the smallest Mol could be related to 
the one-phonon quadrupole structure. The bands with the intermediate values of the Mol could be based on the 
two-phonon quadrupole excitations. 

If the moments of inertia do indeed carry information on the inner structure of the bands, then the numbers of 
excitations with different structure are comparable. This would be in contradiction with the IBM calculation which 
predicts predominantly the octupole two-phonon structure of 0+ excitations (see below). The nature of 0+ excitations, 
derived from calculations in the framework of the QPM as predominately quadrupole, is in contradiction with both 
the IBM calculation and with the above mentioned empirical observation (see below). 

D. Excited states with spins higher than 0"*" 

Other states, mainly with spins 2+, 4+ and 6+ are intensively excited in the (p,t) reaction. The nature of these 
states may only be assumed. Some of these states could belong to the collective bands based on 0^ states. Some of 
the 2+ states could be quadrupolar (one-phonon) vibrational states with correspondence to 0+ excited states, since 
in deformed nuclei every excitation of angular momentum I'^ splits into states distinguished by their K quantum 
numbers ranging from to /. Some 4+ excited states could be hcxadecapolc vibrational excitations and 2+ and 0"*" 
states should correspond to this class of states for the same reason. If this speculation reflects reality then the number 
of 0+ states has to be the largest. However, the observation is in contradiction with this conjecture unless very weak 
0+ excitations are not seen in the (p,t) reaction. Attributing the underlying structure to each of the observed states 
is not possible with the presently available experimental data. To give at least a hint for the structure of these states 
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and corresponding collective bands. 



the experimental energy distribution of the (p,t) transfer strength for the 0"*", 2+ and 4+ excitations is plotted in 
Fig. 5. 

A most remarkable feature in Fig. 5 is the most strongly excited 4+ state at 1337 keV, which can be related to 
the strongly excited 0+ state at 635 keV and 2+ state at 781 keV. The inertial parameters derived from the bands 
based on these states are practically the same. The only explanation is that these states have the same structure. 
One can assume that these states arc a triplet originating from an excitation of multipolarity 4 as in the case of 
a quadrupole two-phonon excitation. However, a corresponding strongly excited one phonon quadrupole excitation 
with spins 0+ and 2+ are not observed at lower energy. Indeed, the 0+ state at 1297 keV, which was identified as 
the /3- vibrational state [18], is rather weak. Thus quadrupole two-phonon excitations have to be excluded. On the 
assumption that these states have a two-phonon octupole structure a quadruplet of states has to be observed with 
spins from 0"*" to 6"*". Experimentally a state with spin 6"*" is identified at 1653 keV, although its excitation is only of 
0.9 /xb, i.e. much weaker than the other states. If this assumption is correct, then the interpretation of the first 0"*" 
state as a two-phonon excitation obtains some confirmation. This quadruplet of states and the corresponding bands 
are displayed in Fig. 10. Note however, that the moments of inertia derived from the band based on the states of the 
multiplet are somewhat smaller than the ones derived for the bands based on negative parity states (Table V). 

An indirect confirmation of the above assumption would be the existence of another quadruplet. The 1~ levels at 
508 keV and 952 keV and the 2~ level at 1079 keV were interpreted as members of a quadruplet of octupole shape 
oscillations with = 0~ to 3~ [18]. The 3~ level at 1259 keV identified from the (p,t) reaction in the present study 
can be the missing member of this quadruplet, as shown in Fig. 10. The energy separation between the K'^ = 0~ 
and JC^ = 1~ and between the JC^ = 1~ and JC^ = 2~ band heads differ strongly. This was explained by a strong 
coupling of the bands based on the last two states [18]. The energy separation between the 3~ level at 1259 keV and 
the 2^ level is close to the one between the second 1~ and 2~ levels, which can also be attributed to the coupling of 
these three states. 
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III. DISCUSSION 
A. 0"*" excitations 

The importance of pairing in the enhancement of the cross section of the 0+ two-nucleon transfer reaction was noted 
already in an earher pubhcation [27, 28]. The superfluid ground states of deformed nuclei are strongly populated in 
the (p,t) reaction due to the large overlap of the wave functions between nuclei with neutron number N and TV ± 2. 
The excited 0+ states can be populated in the (p,t) reaction due to the fluctuation of the pairing field. Such states 
represent pairing vibration modes [27, 28]. It was realized that pairing correlations in deformed nuclei are induced 
not only by monopolc but also by quadrupolc pairing interactions. If both pairing fields are comparable in strength 
and if there is a nonuniform distribution of oblate and prolate single-particle orbitals around the Fermi surface, this 
may give rise to 0+ excitations treated as pairing isomers decoupled from the deformed superfluid ground state [24? 
]. An asymmetry between (p,t) and (t,p) cross sections predicted by this model was confirmed in an experiment by 
Casten et al. [30]. 

The structure of excited 0^ states in deformed cvcn-cvcn nuclei is still a matter of controversial discussion despite 
intensive investigation. Traditionally the first excited Ot state has been interpreted as the beta- vibrational excitation 
of the ground state. However, in many nuclei the state has only weak transitions to the ground-state band, 
while strong electric quadrupole transitions to the gamma band have been found [31]. This contradicts the traditional 
interpretation, since a transition from a beta- vibrational state to the gamma band is suppressed due to the destruction 
of a beta phonon, and, at the same time, the creation of a gamma phonon. The unclear situation led to an intense 
debate about the structure of low-lying O"*" states. 

Maher et al. [2] were the first who noticed an interesting feature of 0^ excited states in the actinides. The strong 
excitations of the first excited 0"*" states in the (p,t) reaction, combined with all other available evidence (rather 
weak E2 transitions to the ground-state band, strong a decays leading to them, the strong Coulomb excitation of 
the associated collective bands) suggest that these states represent a new and stable collective excitation, different in 
character from the most common formulation of the pairing vibration as well as from the /? vibration usually found in 
the deformed rare-earth nuclei. The second excited 0+ states in actinides (firmly assigned) demonstrate completely 
different features [18, 32, 33]. Weak excitation in the (p,t) reaction, relatively strong E2 transitions to the ground- 
state band and a small B{E1)/ B{E2) ratio for transitions to 2+ and 1~ states give evidence that they could be the 
usual (3 vibrational states. For ^'^°Th this is the case for the level at 1297.8 keV [18]. Otsuka and Sugita [34] applied 
the spd/-interacting boson model to the actinide nuclei, aiming to get a unified description of quadrupole-octupole 
collective states. They suggested the first excited O"*" band be referred to as the "super (3 band" thus emphasizing 
the difference of the structure of this state from that of the usual beta-vibrational state. They also predicted the 
existence of the second excited 0+ band to be the usual (3 band that was confirmed later. There are evidences for the 
first 0-|- state in ^^^Th (quadruplet of states, large moment of inertia comparable to that for the octupole- vibrational 
bands, together with the features noticed by Maher et al. [2] and listed above) to carry the two-phonon octupole 
nature. However, the B(E1)/B(E2) ratio for transitions to the 1- and 2+ states is even smaller for this state than for 
the state at 1297 keV: ~ 4 • 10"'' fm~^ compared to ~ 7 • 10"'^ fm~^ [18]. Intuitively one would expect that a large 
B(E1)/B(E2) ratio might be characteristic for a two-octupole-phonon excitation, whereas a small ratio might indicate 
a P shape oscillation. That is true for the state at 1297 keV (which indeed is a beta vibrational state), but not true 
for the first excited state. Moreover, the IBM and the QPM predict one phonon quadrupole nature for this state (see 
below). Therefore the available data do not allow the firm conclusion on the nature of this state. 

Understanding of the structure of the higher excited states remains a challenge for further experimental studies 
(e.g. 7 spectroscopy in the (p,t7)-reaction) and for nuclear theory. The first attempt to explain experimental data 
of a large number of O"*" excited states in ^^^Gd [35] was a phenomenological approach [36] based on the extended 
interacting boson model {spdf-lBM), which accounted for a large fraction of the observed states. The importance of 
the octupole degrees of freedom was revealed. The first microscopic approach was performed in the framework of 
the projected shell model (PSM) [37], using a restricted space spanned by two and four quasiparticle states. The 
IBM calculation reproduced satisfactorily all energy levels in ^^*Gd and gave small E2 decay strengths for them. 
Soloviev and co-workers [38, 39] applied the quasiparticle-phonon model (QPM) to get a microscopic understanding 
of low-lying 0+ states. The QPM was also applied to -'^^^Gd [40]. It predicts an sizable fraction of the 0+ states to have 
large or dominant two-phonon components, mainly built from collective octupole phonon components, in agreement 
with the IBM calculations [36]. 

The resiilts of a specific analysis of experimental data for actinide nuclei were compared to the .spd/-IBM calculations 
in our previous paper [3]. After publication of these data, Lo ludice et al. carried out a calculation for these nuclei 
in the framework of the QPM [41]. Although energies, E2 and EO transition strengths and two-nucleon transfer 
spectroscopic factors were computed, only some of them could be compared with the then available experimental 
data. We compare our new data for 230rpjj ^j^j^ ^j^g results of both calculations below. 
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No extended calculations have yet been carried out for the excitation of the two-quasiparticle (2QP) modes, which 
are expected to occur at excitation energies of about twice the pairing gap energy, i.e. about 1.5 McV in our case. 
Only very restricted microscopic calculations for low energies were attempted by Ragnarsson and Broglia [24]. It 
would be desirable to extend these calculations to higher energies. Just above 1.5 MeV, exactly in the range 2.0-2.5 
MeV, a bump in the distribution of the (p,t) transfer strength is observed (Fig. 5). At least some of these excitations 
could be of 2QP-nature. A model by Rij and Kahana [42] describing the 0+ state as a pair of holes in the oblate 
1/2[501] Nilsson level should also be mentioned. 

The monopole pairing vibration state for neutron-pair excitation (n-MPV) is expected to be strongly excited in the 
(p,t) reaction because of the large overlap of the wave function of such a state with that of the target nucleus ground 
state. A dominant (p,t) cross section for a single state in the spectrum (besides the ground and the first excited 
states) is observed in some nuclei. In ^^^Pa the cross section for the L = transfer to the state at 1500 keV is about 
15% of that for the ground state and comparable to that for the first excited 0+ state in ^"^"Th [4]. However, this 
behavior is unexpected, since this energy corresponds more to that of the proton-pair excitation (p-MPV), where the 
cross section is expected to be much weaker. At the same time, a dominant cross section in ^^*Th is observed for a 
single state at about 2.1 MeV, much higher than in ^^^Pa, which can be considered as being due to the ^^^Th core 
plus a proton. These facts, as well as practically no correspondence of the energy distribution of the (p,t) strength in 
these two nuclei, need a theoretical explanation. No dominant excitation at higher energies is observed in ^'^"Th. In 
the case of a relatively dense spectrum of 0+ states, fragmentation of the n-MPV state to nearby states is possible. 
Such a group of states in "^^^Th around 1600 keV could be a result of such fragmentation. The summed cross section 
of this group is about 13% of that for the ground state. A similarly dense spectrum occurs also in ^^**Th and ^^^Pa; 
nevertheless, these nuclei demonstrate dominant excitation of individual states. Unfortunately, no calculations have 
been carried out for odd nuclei, though they are planned for ^^^Pa [43]. There is no clear picture concerning the 
problem of the MPV as well as of the 2QP states. 
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FIG. 11: Lowest 0^ bands identified earlier and suggested in this analysis in comparison to the spdf-lBM calculations. 



B. IBM calculations 



Although ^'^"Th is considered a vibrational-like nucleus, the inclusion of the octupole degree of freedom in the 
description of its properties turned out to be important [18]. The role of the octupole degree of freedom in deformed 
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actinide nuclei and the related description with / bosons added to the IBM in the sd boson space {sdf-lBM) has 
been studied in [45] . Despite reproducing reasonably well the main features of the observed low- lying negative parity 
states in the rare earth nuclei [44], the sdf-lBM was not so successful for the actinide nuclei. A better reproduction 
of the relevant data is obtained if a p boson is included in addition to the / boson without seeking an understanding 

of its physical nature [46] . 

The IBM Hamiltonian in the spdf space includes a vibrational contribution and a quadrupole interaction in the 
simple form [36, 46]. 

H = edfid + ^pfip + effif - nQspdf ■ Qspdf, (1) 

where , Cp, and e/ are the boson energies and n^, fip, and fif are the boson number operators. Note that the same 
strength k of the quadrupole interaction describes the sd bosons and the pf bosons. The Qspdf quadrupole operator 
has the form 

Qspdf = Qsd + Qpf = [s^d + dh]^'^^ 

- Xsd[d^^^^ + (3/5)\/7[pt/ + /tp](2) 

- (9/10)x/3[ptp](2) - (3/10)x/42[/t/](2) (2) 

The factor in front of the [rf^rf]^^^ term may be adjusted by introducing an additional parameter Xsd- 

The IBM parameters in the sd boson space are determined by the low energy spacing of the ground state band 
and the /'^ = 2^,K'^ = 0+ and I'^ = 2^,K'^ = 2+ band heads, respectively. The pf boson parameters are chosen 
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to reproduce the K'^ = 0~ and K"^ = 1~ band heads; they are determined by the experimental energies of the 
I'" = 1^,K^ = 0~ and I'^ — l^jK^ = 1~ states. The total number of bosons is 11. The number of negative parity 
bosons is allowed to range from to 3. Relative to the calculations presented in our earlier paper [3], only some 
parameters are slightly changed in the present calculations: ep=-1.00, ed=0.25, e/=-0.90, k=-0.014, Xs<i=-l-00. As 
for rare earth nuclei [36], in the spdf-lBA calculations mixing between d and pf bosons is neglected and the / (and 
p) bosons account for octupole collectivity. In Fig. 11 we display the firmly assigned experimental 0+ states and the 
sequences on top of them, in comparison with the spdf-lBM calculations. 

The full experimental spectrum of the 0+ states in ^^°Th, including relatively firm and tentative assignments, and 
the results of spdf-lBA calculations arc compared in Fig. 12. In the energy ranges covered experimentally, the IBM 
predicts seven excited 0+ states of pure sd (quadrupolar) bosonic structure, and twelve excited 0+ states which have 
two bosons in the pf boson space. They could be related to octupole two phonon excitations (OTP). The nature 
of the states according to the calculations is indicated in Fig. 12. The number and the nature of the 0+ states in 
the IBM calculations depends on the parameter values but the solution presented in Fig. 12 is stable for a realistic 
variation of these parameters. A complete study of this dependencies would be instructive but it is not the goal of 
the present paper. 

We can see a reasonable correlation in excitation energy between experiment and calculation up to 1.4 MeV. The 
calculation predicts the first excited 0+ state as a quadrupole excitation in the sd space at an energy close to the 
experimental state at 635 keV. The second excited 0+ state is predicted as an octupole two-phonon excitation and 
could correspond to the tentative assigned 0+ state at 953 keV. The following two are predicted to have a quadrupole 
structure and correspond energetically to the tentative assigned 0+ state at 1126 keV and to the 0+ state at 1297 
keV, respectively. The latter was, indeed, considered the /? -vibrational state in Ref. [18]. 

Thus the IBM calculation with the parameterizations used predicts for "^^^Th 20 excited 0+ states in the energy 
range below 2.7 MeV. Accounting in addition for the presence of monopole pairing vibrational states, two-quasiparticle 
states and perhaps a state from hexadecapole collectivity, not included in the calculation, we can consider 24 observed 
0"*" excitations as nearly perfect agreement with calculated number of such excitations. But there is no clarity 
concerning the nature of these excitations without additional experimental information. The IBM also fails completely 
to reproduce the (p,t) spectroscopic factors. The calculated first excited 0+ state occurs with about one percent of 
the transfer strength of the ground state, and the higher states are even weaker, whereas experimentally the excited 
states show at least about 70 % of the ground state transfer strength. 

A calculation in the framework of the spdf-lBM gives 20 levels of spins 2+, 4+ and 6+ in comparison with 40, 32 
and 11 identified levels of these spins in the energy region below 2.7, 2.3 and 2.4 MeV, respectively. The spectrum of 
2+ and 4+ states in comparison with the spdf-lBM calculation is given in Fig. 12. Most of the 6+ states could not 
be observed in the (p,t) reaction because of their low cross section. As far as 2+ and 4+ states are concerned, the 
number of experimental levels is much higher than the prediction of the spdf IBM model. 



The ability of the QPM to describe multiple 0+ states (energies, E2 and EO strengths, two-nucleon spectroscopic 
factors) was demonstrated for ^^*Gd [40]. An extension of the QPM to describe the 0+ states in the actinides [41] 
was made after our publication on the results of a preliminary analysis of the experimental data [3]. The present 
calculations aim to explain the results of the detailed analysis of the experimental data for ^^'^Th. 

In the QPM [47] the Hamiltonian in a separable 

generalized form is adopted to generate the quasiparticle RPA phonons described by the operators 



9192 

The Hamiltonian expressed in terms of these phonon operators is diagonalized in the space spanned by one- and 
two-phonon states. The QPM eigenstates have the structure 



where XK label the multipolarity and magnetic component of the phonon operator. Each of these states represents 
the intrinsic component of the total wave function 



C. QPM calculations 
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FIG. 13: Comparison of experimental and calculated (QPM) 0"'", 2+ and 4+ relative level reaction strengths for the (p,t) 
reaction. The 2+ and 4+ states assumed to belong to rotational bands (see Sec. II C) are not included. The experimental 
strengths for Og.g., 2^ and 4^ being rotational states are normalized to 10 (the latter two are not shown in figure). 



where i?£f^ is the Wigner matrix. No free parameters are used in these calculations, all physical input quantities 
and constants are determined by an independent fit to the experimental data in neighboring odd nuclei. For details 

see [41]. 

The experimental spectra of the 0+, 2+ and 4+ states in ^^°Th are compared with the results of the QPM calculations 
in Fig. 13. The QPM considers only vibrational 2+ and 4+ excitations in the even-even nuclei, therefore all other 
excitations have to be excluded from the comparison. Only the 2+ and 4+ states not belonging to rotational bands 
(see Sec. II C) are included in Fig. 13, i.e. these states are assumed to be mainly of vibrational structure. The QPM 
generates 23 0+ states below 2.8 MeV in fair agreement with the 24 identified states. At the same time the QPM yields 
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less 2+ and 4+ excited states than the ones observed experimentally. The numbers of 2+ and 4+ states generated 
up to 3 MeV are 20 and 17, respectively, considerably less than the observed 50 and 27 corresponding experimental 
levels (Fig. 13). It seems that taking into account the two-quasiparticle and pairing vibrational states not excluded 
from consideration (there are problems in their identification) cannot explain this large difference. Beside the ground 
state, experiment reveals one very strong peak for every of G+, 2+ and 4+ state and small strengths for other states. 
The calculation correctly yields strong peaks close in magnitude and position to the experimental ones for the 0+, 2+ 
and 4+ states. These strong peaks form the suggested multiplet shown in Fig. 10. The calculation correctly yields 
small strengths for other peaks in the case of 0"*" excitations. At the same time, besides the first strong peaks, the 
QPM predicts for the 2+ and 4+ states two other strong peaks not observed experimentally. 

In contrast to the spdf-lBM, the QPM is able to reproduce the two-neutron transfer strength. The wave functions 
(5) can be used to compute the (p,t) normalized transfer spectroscopic factors 



Sn{P,t) 

where the amplitudes are given by [48] 



rn(p,t) 



(6) 



r„(p,i) = {^l^j^j^N -2\J2r'YiKag,ag,\^o,N). (7) 

The amplitude To{p,t) refers to the transitions to the / members of the ground state rotational band [41]. In Fig. 14 
we present the increments of the (p,t) strength to the 0"^, 2+ and 4+ states and that of the spectroscopic factors 
derived from the DWBA analysis for these states. They are given relative to those for corresponding states of the 
ground state band and are compared with the calculated normalized spectroscopic factors. Since the DWBA analysis 
for the 0''' states included different configurations of the transferred neutrons, only the (p,t) strength ratio is given 
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in Fig. 14 for these states. For other states only the (299/2)^ neutron configuration was accepted. According to the 
DWBA formalism in the case of the same neutron configuration and the direct one-step transfer the strength ratios 
have to be close to the spectroscopic factor ratios. This is not the case for ^^°Th, since a considerable difference of 
these two ratios is observed. 

As one can see, the calculations of 0+ for the (p,t) strength ratio are in fair agreement with the experiment. For the 
2+ states the calculations are in good agreement with the measured spectroscopic factor ratios, however almost two 
times smaller than the measured strength ratios. We have to note that the (p,t) strengths for the 0+ and 2+ states not 
firmly assigned are small and do not influence considerably the results of comparison . At the same time withdrawal of 
these states from consideration will only improve the agreement with the calculations. For the 4"*" states the calculated 
spectroscopic factors are more than two times larger than the experimental strength and spectroscopic factor ratios. 
The difference possibly stems from the change of the (p,t) cross section caused by the inclusion of additional two-step 
paths of the neutron transfer for the 2+ and 4+ members of the ground state band. As one can see from Fig. 7, 
the angular distributions for these states differ considerably from those for direct one-step transfer and can be fitted 
only by an inclusion of a two-step excitation number which is the largest for the 4+ state. For the 0+ states, where 
only the one-step transfer is possible, the agreement between calculation and experiment is good. For both 2+ and 
4+ states the energies at which the experimental strength first appears agrees with the calculated ones. This means 
that the energies of the first vibrational 2+ and excitations predicted by the QPM are in fair agreement with the 
experiment (recall that the states belonging to the rotational bands are excluded from this comparison). 

The nature of 0+ excitations as well as for 2+ and 4+ states in the QPM differs greatly from that in the s;>ci!/-IBM. 
In all low-lying states quadrupole phonons are dominant and the octupole phonons arc predicted to play a relatively 
modest role. This might be acceptable in the vibration-like ^^°Th, but the same is predicted for the octupole soft 
^^^Th. The spectrum is explained from the QPM calculation procedure (Pauli principle) as a redistribution of the 
strength of the lowest two-octupole phonons among many closely packed QPM 0+ states [41]. To assess this aspect, 
calculations for octupole deformed lighter isotopes of Th would be important. 

Other predictions of the QPM can be tested experimentally only for the lowest states. The measured values of 
B{E2;0f — » 2g ) = 1.1 W.u. is close to the computed value of 1.7 W.u. At the same time the calculated monopole 
transition strength p'^{EO;0'^ 0+) = 1.48 • 10~^ is two orders of magnitude smaller than the experimental value 
of 126(13) • 10^'^ [18]. Again, we would like to stress the necessity of systematic measurements of electromagnetic 
properties of higher excited 0"*" states to understand their nature. 

IV. CONCLUSION 

Excited states in ^^°Th have been studied in (p,t) transfer reactions. About 200 levels were assigned using a 
DWBA fit procedure. Among them, 24 excited 0+ states have been found in this nucleus, most of them have not 
been experimentally observed previously. Their accumulated strength makes up more than 70 % of the ground state 
strength. Firm assignments have been made for most of the 2"*" and 4'*" states and for some of the 6^ states. These 
assignments allowed the identification of sequences of states which have the features of rotational bands with definite 
inertial parameters. The 2+, 4+ and 6+ states not included in these bands have been considered as vibration-like 
excitations. Multiplets of states arc suggested which can be treated as one- and two-phonon octupole quadruplets. 
The experimental data are compared with spdf-lBM and QPM calculations. Giving an approximately correct number 
of 0+ states, these models provide different predictions for the structure of these states. They are also in conflict 
with the apparent structure of the states inferred from the moments of inertia of the rotational bands built on them. 
More specifically, as follows from the moments of inertia, the 0+ states have a different intrinsic structures, which is 
in contradiction with the predictions of both the IBM (predominantly octupole bands) and the QPM (predominantly 
quadrupole bands). A remarkable feature of the QPM is the prediction of strong first vibrational excitations close 
in magnitude and position to the experimental ones. The numbers of 2+ and 4+ states are underestimated by both 
theories. Spectroscopic factors from the (p,t) reaction, and the trend in their change with the excitation energy, 
are approximately reproduced by the QPM for the O"*" and 2+ states and overestimated by theory for the 4+ states. 
The lack of additional information does not allow for final conclusions on the validity of the theoretical approaches. 
Therefore we hope that our new data will stimulate further experimental and theoretical studies. Accurate experiments 
and a detailed analysis similar to the present work are desirable for other nuclei in this region. Challenging experiments 
on gamma spectroscopy following (p,t) reactions would give much needed information. 
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APPENDIX: Angular distributions and their fit by the CHUCKS 
calculations for all excited states observed in the present study 



8 
6 
4 
2 


3 
2 
1 


4 
3 
2 
1 

0,8 
0,6 
0,4 
0,2 


3 
2 
1 


2 
1 


1,5 
1 

0,5 


2 

1 





5 1762.3 keV 
^5 




_i 1 i_ 




3 F 

2 - 

1 - 







' 1 ' 1 ' 1 

1851 .4 keV 






1 


I.I.I 




1859.3 keV 



' 1 ' 
1 


' 1 ' 1 
1939.8 keV " 

1 1 




_ ' 

- 2 ^ 
1 


■ 1 ■ 1 _ 

1947.0 keV . 
1 1 



10 20 30 40 



10 20 30 40 

angle 



10 20 30 40 



12 
9 
6 
3 


4 

2 


1 


0,8 
0,4 



2 

1 



8 
6 
4 
2 

3 

2 

1 


8 

6 

4 

2 







15 
10 
5 


2 

1 


2 
1 




"T ' 1 ' 1 ' [— 

2226.0 keV 




1 ' 1 ' — r- 

2241 .0 keV 




_i 1 I 



1 ' r 

2249.9 keV 




10 20 30 40 




20 
15 
10 
5 


4 





10 



2305.4 keV 




— ' 1 ' 1 ' 1 ' 1 : 

2311 .2 keV - 



■ s 


2317.7keV 




5 ■ 



5 2329.6 keV 
"5 




10 20 30 40 

angle 




i 1 i 1 i 




' 2383.8 keV 

I.I. 






10 20 30 40 



4 

3 

2 

1 


6 

4 

2 


4 

3 
2 
1 


10 

5 


4 

2 



8 
6 
4 
2 

8 
6 
4 
2 

3 

2 

1 









2824.4 keV 







10 20 30 40 

angle 



2834.0 keV 
5 







2870.6 keV 







10 20 30 40 



6 
4 
2 



3 

2 

1 



8 

6 

4 

2 



10 
8 
6 
4 
2 


6 
4 
2 


2 
1 



4 
3 
2 
1 

10 
8 
6 
4 
2 




"T ' 1 ' 1 ' T" 

2906.4 keV 




"T ' 1 ' 1 ' r 

3009.9 keV 



2913.6 keV - 





2950.5 keV 




ii 




10 20 30 40 





3072.6 keV 



"T ' 1 ' 1 ' [— 

3083.8 keV 




10 20 30 40 

angle 




3135.9 keV 



3147.4 keV 

I 




10 20 30 40 



